#include<iostream>
#include<fstream>
#include<cmath>
#include<vector>
#include "../programming1/EquationSolver.h"
#include "../HW3/Thomas.h"




 int main(){
    vector<double> d={2,2,2};
    vector<double> l={1,1};
    vector<double>  u={-1,-1};
    vector<double>  b={1,2,3};
    // int n = d.size();
    vector<double>  M(3);
    M=Thomas(d,u,l,b);
    for (int i=0;i<M.size();i++)
        cout<<M.at(i)<<endl;
    return 0;
 }